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High-precision data from observation of the cosmic microwave background and the large scale 
structure of the universe provide very tight constraints on the effective parameters that describe 
cosmological inflation. Indeed, within a constrained class of ACDM models, the simple \4>^ chaotic 
inflation model already appears to be ruled out by cosmological data. In this paper, we compute 
constraints on inflationary parameters within a more general framework that includes other phys- 
ically motivated parameters such as a nonzero neutrino mass. We flnd that a strong degeneracy 
between the tensor-to-scalar ratio r and the neutrino mass prevents A(/>'* from being excluded by 
present data. Reversing the argument, if \(f>* is the correct model of inflation, it predicts a sum of 
neutrino masses at 0.3 — » 0.5 eV, a range compatible with present experimental limits and within 
the reach of the next generation of neutrino mass measurements. We also discuss the associated 
constraints on the dark matter density, the dark energy equation of state, and spatial curvature, 
and show that the allowed regions are significantly altered. Importantly, we find an allowed range 
of 0.094 < Q.ch^ < 0.136 for the dark matter density, a factor of two larger than that reported in 
previous studies. This expanded parameter space may have implications for constraints on SUSY 
dark matter models. 



I. INTRODUCTION 

The past few years have seen a dramatic increase in 
the precision of cosmological data, ranging from mea- 
surements of the cosmic microwave background (CMB) 
anisotropics by the Wilkinson Microwave Anisotropy 
Probe (WMAP) satellite [J, i, i] and the large scale 
structure (LSS) of the universe by the Sloan Digital Sky 
Survey (SDSS) 0, to the observation of distant type 
la supernovae (SNIa) @. All these measurements point 
to the so-called concordance model of cosmology, wherein 
the physical parameters are the baryon density Qt, the 
matter density 17™, the dark energy density and the 
present Hubble expansion rate Hq. The model geometry 
is flat so that Q/^ = 1 — flm, and the initial perturbations 
are assumed to be adiabatic and Gaussian, with a power 
law spectrum described by a spectral index Us and an 
amplitude As. Together with the optical depth parame- 
ter T, this six-parameter "vanilla" model provides a good 
fit to all observational data to date. 

A common assumption in cosmological parameter es- 
timation is that one can always improve a fit marginally 
by including extra free parameters. This assumption has 
led to the adoption by many authors of the Occam's ra- 
zor approach, in which an extra parameter is retained 
only if by its inclusion the goodness-of-fit of the model is 
substantially improved. Indeed, the success of the vanilla 
model is rooted in the fact that, given the current data, 
no addition of a single extra parameter produces a 
value that is significantly lower. 

However, there are many more physically well mo- 
tivated parameters beyond the vanilla model. Indeed, 
some of these, such as a nonzero neutrino mass, are 
known to be present. In such cases, a blind enforcement 
of Occam's rule can lead to significant underestimation 



of parameter errors, as well as bias in the parameter es- 
timates. One well-known example is the interplay be- 
tween the dark energy equation of state and the neu- 
trino mass 0, B [1| • When the dark energy equation of 
state is allowed to vary, the neutrino mass bound is re- 
laxed by almost a factor of three if only CMB and LSS 
power spectrum information is used. Conversely, by im- 
posing a prior on the neutrino masses accordingto the 
Heidelberg-Moscow claims [HI Ell ) reference finds 
that a cosmological constant is ruled out at more than 
95 % C.L. by CMB-fLSS-l-SNIa data. 

One could argue that parameter estimation coupled 
with Occam's rule is a "bottom-up" approach, for which 
a full Bayesian analysis complete with Bayes factor calcu- 
lations may also be appropriate However, if one's 
aim is to exclude specific models, then a more conserva- 
tive approach that takes into account possible degenera- 
cies between the "standard" and the "new" parameters is 
warranted. Such a "top-down" approach does not neces- 
sarily imply a decrease in the predictability of the model. 
In fact, we will show that given the present cosmological 
data, a nonvanishing neutrino mass could be viewed as a 
prediction of the A(/)^ inflationary model. We argue that 
when constraining or excluding specific theoretical mod- 
els, one should in principle allow for uncertainties in all 
physically well-motivated parameters, even if they have 
a priori no direct link to the models concerned. If, for 
instance, it turns out later that the universe is indeed 
composed of a nonvanishing neutrino fraction, it would 
be counterproductive to have already discarded a model 
of inflation that predicts this outcome. 

In the present work, we investigate in this spirit how 
parameter constraints change when the parameter esti- 
mation analysis is performed within a much more general 
model framework. In principle there are some twenty 
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or more parameters that could influence cosmology, al- 
though the precision of the present data is not yet suffi- 
cient to probe some of them (e.g., the primordial helium 
fraction and the effective sound speed of dark energy). 
Here, we focus on a 11-paramctcr model outlined below. 



A. The model 

We test a general 11-parameter model space consisting 

of 

6 = {ujc.ujb, U,^k,w,Ho,ns,r,as,T,As). (1) 

The vanilla model is defined by fi, = Qk = f = <^s = 0, 
and w = — 1. In addition, we marginalise over a nuisance 
parameter b which describes the relative bias between the 
observed galaxy power spectrum Pg{k) and the underly- 
ing dark matter spectrum Pc(fc) via Pc{k) — lP'Pg{k). 

Three different parameter sets will be considered in 
this work: 

• Set A: All 11 parameters. 

• Set B: A 10-parameter set with fifc = 0. This pa- 
rameter set covers all standard inflationary models 

• Set C: A 9-parameter set with 0,^ — cts — 0. This 
reduced set corresponds to the large subset of the 
zoo of inflationary models that predict negligible 
running, including large field chaotic inflation mod- 
els 



1. Matter content 

We assume the matter content to be specified by the 
following parameters: the curvature fifc = 1 — ^l„i — JIa, 
the physical dark matter density = ^ch^, the baryon 
density LUt = Ob/i^, the neutrino fraction = fl^/ftc, 
and the dark energy equation of state parameter w. 

Other parameters not included here, but which could 
have an observable effect, include a time-dependent dark 
energy equation of state, nonstandard interactions in any 
of the dark sectors (cold dark matter, neutrinos, or dark 
energy), etc. We mention this as a caution that while our 
parameter space is much larger than that normally used 
in parameter estimation analyses, it is not necessarily 
complete. 
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FIG. 1: Two-dimensional 68 % and 95 % C.L. contours for the 
inflationary parameters n^, r, and as, using the full data set 
and parameter set B, and marginalised over the other (10 — 2) 
parameters. 



2. Initial conditions 

The initial conditions for structure formation are as- 
sumed to be set by inflation, characterised by scalar and 
tensor fluctuations with amplitude As and At = rAs re- 
spectively. Each component is specified by a spectral 
index or rit, and the inflationary consistency relation 



requires that rit ^ —r/8. However, the precision of cur- 
rent data is not yet at a level where a violation of the 
consistency relation can be tested. This also means that 
while the running parameter should be included for 
the scalar spectrum, the inclusion of its tensor counter- 
part at would have no effect. This set of initial param- 
eters encompasses all standard inflationary models, but 
not models with features from potential steps, particle 
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FIG. 2: Degeneracies between r and fv,w for the full data set 
and parameter set B, marginalised over (10 — 2) parameters. 



production, etc. during inflation. We deflne as at the 
pivot scale k = 0.002 Mpc^^, in concordance with most 
recent analyses. 

Note that an alternative approach would be to perform 
the analysis directly in terms of the slow-roll parameters 
instead of the observables , as, and r [l^ [H, [13] ■ Par- 
ticularly for models where Ug is not negligible this can 
lead to somewhat different results. However, in models 
with small as, such as chaotic inflation, the results are 
identical. 



TABLE I: The 95 % C.L. allowed ranges for n^, r and aa for 
parameter set B, marginalised over the other (10 — 1) param- 
eters. 



Parameter 


WMAP-fSDSS 


Full data set 


Us 


0.97 ^ 1.35 


0.98 ^ 1.28 


r 


^ 1.05 


^ 0.81 


Qs 


-0.140 ^ -0.005 


-0.135 ^ -0.004 



isation over the bias b and the nonlinear correction pa- 
rameter Qnl- 

Baryon acoustic oscillations (BAO) In addition to 
the power spectrum data we use the measurement of 
baryon acoustic oscillations in the two-point correlation 
function The analysis is performed following the 

procedure described in [13, Hlj (see also including 
analytic marginalisation over the bias 6, and nonlinear 
corrections with the HALOFIT ^] package. 

Type la supernovae (SNIa) We use the luminosity 
distance measurements of distant type la supernovae pro- 
vided by the Supernova Legacy Survey (SNLS) Q . 

Lyman-a forest We do not include data from the 
Lyman-a forest in our analysis. These data were used 
in some previous studies that found very strong bounds 
on various cosmological parameters 



2J|. However, the 



strength of these bounds is due mainly to the fact that the 
Lyman-a analysis used in (24| leads to a much higher nor- 
malisation of the small-scale power spectrum than that 
obtained from the WMAP data. Other analyses of the 
same SDSS Lyman-a data find a lower normalisation, 
in better agreement with the WMAP result [H, [H, [l^ • 
This kind of discrepancy between different analyses of the 
same data probably points to unresolved systematic is- 
sues, and for this reason we prefer to discard the Lyman- 
a data entirely. 

For a large part of our analysis we use two different 
combinations of data sets, one consisting of WMAP and 
SDSS data only, and one which uses in addition data from 
SNIa (SNLS) and BAO. The latter case is sometimes 
referred to as "the full data set" . 

We perform the data anal ysis using the publicly avail- 
able CosmoMC package [28l l29l|. modified to include the 
BAO likelihood calculations. 



II. DATA ANALYSIS 

Cosmic microwave background (CMB) We use CMB 
data from the WMAP experiment after three years of 
observation [l|, |3, [!]■ The data analysis is performed 
using the likelihood calculation package provided by the 
WMAP team on the LAMBDA homepage [l9|. 

Large scale structure (LSS) The large scale structure 
power spectrum of luminous red galaxies (LRG) has been 
measured by the Sloan Digital Sky Survey (SDSS). We 
use the same analysis technique on this data set as advo- 
cated by the SDSS team 0,11], with analytic marginal- 



Ill. RESULTS 
A. Inflationary parameters 

Almost all inflationary models predict fifc to be zero. 
This prediction is also supported by our analysis of pa- 
rameter set A (see Sec. IIIIDj) . Therefore, in this section, 
we will work with the reduced 10-parameter set B, in 
which rife is already fixed at zero. Figure [1] shows the 2D 
likelihood contours for the parameters Us, r and as using 
the full data set and parameter set B. These contours are 
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FIG. 3: Two-dimensional 68 % and 95 % C.L. contours for 
and r, using parameter set C (consistent with predictions of 
chaotic inflation), and marginalised over (9 — 2) parameters. 
The upper panel uses WMAP+SDSS data and the lower the 
full data set. The two short black/solid lines with boxes at 
the ends correspond to predictions of A<^'* (top left) and m^0^ 
models of inflation, with 46 to 60 e- foldings (left to right). 



obtained by marginalising over the other (10 
eters not shown in the plot. 



param- 



Figure [T] should be compared with, e.g., Figs. 2 and 3 
of Kinney et al. [13], which use data from WMAP and 
SDSS, and a parameter set similar to our set B but with 
fu and w fixed at and —1 respectively. The compar- 
ison reveals that the two sets of likelihood contours are 
roughly similar, but with one important exception: the 
allowed range for the tensor-to-scalar ratio r in our case 
is much larger even in the light of additional data! 

In order to understand this effect we plot in Fig. [2] the 
2D likelihood contours for r and our additional parame- 
ters fn and w. Interestingly, a substantial degeneracy ex- 
ists between r and the neutrino fraction which in turn 
allows r to extend to much higher values. Table [J displays 
the ID 95 % C.L. allowed ranges for rig, r, and as, as- 
suming parameter set B and using both WMAP-I-SDSS 
only and the full data set. 



B. Chaotic inflation 

Single field inflation models with polynomial potentials 
generally predict negligible running. These models are 
thus represented by our 9-parameter set C in which as = 
0. The corresponding 2D likelihood contours for Ug and 
r, marginalised over the other (9 — 2) parameters, are 
shown in Fig. [3] 

Figure El should be compared with Fig. 4 of Kinney et 
al. [301 , with Fig. 14 of Spergel et al. Jj] , and with Fig. 19 
of Tegmark et al. H (see also [U). In all cases our 
WMAP -f SDSS contours encompass a markedly larger re- 
gion. In particular, even with the inclusion of SNIa and 
BAO data, we find that the simplest At/)'' model is still al- 
lowed by data, contrary to the conclusions of [ll. 13. Isol. [Slj . 
We note that the endpoints of the model lines in Fig. |3] 
correspond to 46 and 60 e-foldings respectively."' Inter- 
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FIG. 4: Degeneracies between between r and ,f^,w for the 
full data set and parameter set C, marginalised over (9 — 2) 
parameters. 



^ When taking into account one-loop effects in the chaotic infla- 
tionary scenario, the model lines in this plot may actually be 
smeared at the percent level [s^l- 
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TABLE II: The ID marginalised 95 % C.L. allowed ranges for Us and r for parameter set C and its subsets. 



Parameter set 


Data set 


Us 






r 


C 


WMAP+SDSS 


0.927 


1.038 


- 


0.51 


C 


WMAP+SDSS+SNLS+BAO 


0.932 ^ 


1.018 


- 


0.41 


c 


WMAP+SDSSlin+SNLS+BAO 


0.931 ^ 


1.025 





0.47 


C, w fixed 


WMAP+SDSS+SNLS+BAO 


0.933 ^ 


1.019 


- 


0.40 


C, w, fi, fixed 


WMAP+SDSS 


0.931 -> 


1.011 


- 


0.31 


C, w, fi, fixed 


WMAP+SDSS+SNLS+BAO 


0.931 ^ 


1.010 


- 


0.30 
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FIG. 5: Two-dimensional marginalised constraints on and 
r for the parameter set C, but with the restriction w = —1. 



estingly, Xcf)* is compatible with data only if the number 
of e-foldings is relatively lar ge, o r equivalently, if the re- 
heating temperature is high [33l Isl]. 

Again, the explanation for our enlarged Hs , r allowed 
region lies in our expanded model parameter space. In 
Fig. m we see that the degeneracy between r and en- 
countered earlier in parameter set B is present also in 
parameter set C, albeit to a smaller extent. If a neu- 
trino fraction of 0.03 — * 0.05 is allowed (corresponding 
roughly to 0-3 ^ 0.5 eV), new parameter space 

opens up for Hs and r. We note in passing that the con- 
verse is not true. Allowing r to run does not change the 
upper bound on the neutrino mass significantly. Inter- 
estingly, this fi, , r degeneracy also means that in its 
simplest form predicts quasi-degenerate neutrino masses 
with a sum in the 0.3 — > 0.5 eV range. This range is com- 
patible with present laboratory limits from tritium beta 
decay experiments, nii, < 2.2 eV [H, [1^, as well as the 
claimed detection of neutrinoless double beta decay, and 
hence detection of the effective electron neutrino mass 
^ee = J2j ^ej'^'^i 0.1 0.9 cV , by the Heidelberg- 
Moscow experiment [13, [HI, ■ The upcoming tritium 
beta decay experiment KATRIN will also probe neutrino 
masses to a comparable level of precision [37]. 

Also of interest is the case of a fixed dark energy equa- 



tion of state w. Figure[5]shows the equivalent of the lower 
panel of Fig. [3] (full data set and parameter set C), but 
with the additional restriction w = —\. Clearly, there 
is very little difference between Figs. [3] and El since the 
combination of SNIa and BAO data effectively fixes w to 
— 1 in the former shown in Fig. [4] 



As a consistency check we present in Fig. [S] also the 
2D constraints on ris , r for the vanilla model with one 
extra parameter r, i.e., the same model analysed in [l|, 
[1, [s^l • The general shapes of the contours in this figure 
are almost identical to those in Fig. 19 in Tegmark et al. 
0] which uses the same data sets. In addition, we find a 
ID 95 % C.L. upper bound of r < 0.31, while Tegmark 
et al. report an almost identical r < 0.33. Kinney et al. 
also found r < 0.31 for the same vanilla+r model [soj . 
but from a combination of WMAP and the SDSS main 
galaxy samples (as opposed to SDSS LRG used in this 
work and in For comparison [sH] found r < 0.26 for 
an analysis of WMAP and 2dF. 



Using additional data from the Lyman-a forest, Sel- 
jak et al. [2^ derived an even stronger upper bound, 
r < 0.22, for the same model space. The reason for the 
improvement is a degeneracy between r and dg, such that 
a higher value of r leads to a smaller preferred value of 
(78 ■ Since the Lyman-a data used in [24!| prefer a high 
value of (jg, a small r value is correspondingly favoured. 
In fact, from a parameter fitting point of view, a nega- 
tive r would be even better. All these conspire to give 
a much stronger upper bound on r. However, as noted 
in Sec. [Hi this phenomenon likely points to a systematic 
uncertainty in the Lyman-a normalisation, rather than a 
genuinely strong constraint on r. 



Finally, we stress again that the difference between 
the allowed rts,r regions in Figs. [3] and [S] lies in a de- 
generacy between r and the neutrino fraction f^,. It 
should also be noted that the addition of SNIa and BAO 
data has very little impact on the vanilla+r model, be- 
cause no strong parameter degeneracies are present in the 
WMAP+SDSS data. With SNIa and BAO included we 
find a ID 95 % C.L. bound of r < 0.30, instead of 0.31 
for WMAP+SDSS alone. 
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C. The effect of non-linearity 



So far we have used exactly the same analysis tech- 
nique as the SDSS team when treating the LRG data. 
However, beyond a wavenumber of approximately k ~ 
0.06 0.07 h Mpc~^, nonlinear effects begin to dom- 
inate the matter spectrum (see, for instance. Fig. 9 of 
[3]). To test whether or not our results are subject to 
these effects, we perform the same analysis as in Fig. [31 
but retain data only up to fc ~ 0.06 h Mpc~^ (band 11). 
We call this reduced data set SDSSlin, and the result is 
shown in Fig. [T] Using only the linear part of the power 
spectrum data has no bearing on our conclusions. In 
fact, the 2D allowed region in ns,r for parameter set C 
is only affected in the region where Us > 1. The SDSS 
data probes Us more precisely when all data points are 
included, and this in turn leads to a truncation of the 
allowed region at high Us- 
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FIG. 7: Two-dimensional marginalised constraints on and 
r for parameter set C. Only the linear part of the SDSS power 
spectrum has been used. 



Table |TI] summarises the ID marginalised constraints 
on Hs and r for parameter set C and its subsets. 



D. Dark matter and dark energy 
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FIG. 6: Two-dimensional marginalised contours for the 
vanilla-l-r parameter space, for both WMAP-I-SDSS only (up- 
per panel) and the full data set (lower panel). 



In order to derive robust bounds on the physical dark 
matter and dark energy properties, all other plausible 
parameters should be allowed to vary. With respect 
to the initial conditions this is almost impossible since 
the most general inflationary models do not necessar- 
ily give smooth, power law-like spectra. Instead, the 
primordial power spectrum can have various features 
BHElHSIillil (see also [M [13, IH, gl, H HH 
for more observationally oriented discussions) , which may 
bias estimates of parameters and their errors. Here we 
present just a small step towards dark matter and dark 
energy parameter estimation in the context of more gen- 
eral models. 

The physical dark matter density flch^ is a crucial 
input in dark matter model building. A prime ex- 
ample of this is models with low energy SUSY where 
the dark matter particle is usually either the neutralino 
or the gravitino. Large regions in parameter space in 
these models have been excluded by the fact that the 
predicted dark matter density is too high or too low 

[IS, il, H [IS ii, y, H, [M!, [si . 

In the vanilla model ^Ich^ is a very well constrained 
quantity, with WMAP-hSDSS giving a 68 % C.L. limit 
of n^h^ = 0.1050t!^;J5|^4j [i]. This corresponds to a rel- 
ative uncertainty of a{nch'^) /VLch"^ ~ 0.04. The SDSS 
collaboration also provide bounds on Vlch? in extended 
models in which one additional parameter is added to 
the vanilla parameter set Q. In most cases the bound 
on Q.ch? does not change significantly. However, when ei- 
ther Vlk or w is allowed to vary, a(Jlch^) /Vtch? increases 
to about 0.06 [4]. 

We have taken this investigation further by calculat- 
ing the bound on Vl^h? for our various parameter and 
data sets. In Fig. [5] we show the joint 2D marginalised 
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constraints on flch'^ and w for three different cases using 
parameter sets B and C. If only WMAP and SDSS data 
are used, a very strong degeneracy between flch^ and w 
weakens the bounds on both parameters. This degener- 
acy is broken when SNIa and BAO are included (as is 
also the case with the degeneracy between and w), 
yielding strong constraints on both parameters. 

When spatial curvature is also allowed to vary, the 
bound on flch^ does change considerably. Figure [H] shows 
the 2D marginalised contours for ilch^^flk and w, us- 
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FIG. 8: Two-dimensional marginalised 68 % and 95 % C.L. 
contours for flch^ and w, using various parameter and data 
sets. Top: Parameter set B, WMAP-I-SDSS. Middle: Param- 
eter set B, the full data set. Bottom: Parameter set B with 
Us — (i.e., parameter set C), the full data set. 



ing parameter set A and the full data set. Here we find 
a{nch?)/nch? ~ 0.1, so that 0.094 < f^^/i^ < 0.136, 
-0.022 < f2fc < 0.026, and -1.19 <w< -0.88 (ID at 
95 % C.L.). It is interesting to compare our more general 
constraints on 51^ with that given by the SDSS collab- 
oration from a vanilla-|-rifc fit (—0.015 < Vtk < 0.023, 
95 % C.L.) [1]; our allowed range is slightly larger even 
in the light of additional data from the distance measure- 
ments of SNIa and BAO. We note also that even though 
the allowed range for Q,ch? increases considerable with 
the inclusion of fifc, the same is not true for the dark 
energy equation of state parameter w. In Table IIIIl we 
summarise the ID 95 % constraints on ^Jt? and w from 
Figs. [5] and El 

Finally, let us we stress that some caution should be 
applied whenever the dark matter density is used as an 
input to constrain models such as the MSSM. Param- 
eter regions that are excluded in the simplest vanilla 
model can easily be allowed in more general models, even 
without the introduction of more exotic features such as 
isocurvature modes. If one is to take one single num- 
ber inferred from cosmological observations as an input 
to constrain particle physics models, then the safest ap- 
proach is to allow for the possibility that cosmology is 
not described by the vanilla model, but by something 
more general. From our calculations, we recommend us- 
ing 0.094 < n^K^ < 0.136 (95 % C.L.), but we caution 
that even this may not be the most conservative estimate. 



IV. DISCUSSION 

We have performed a detailed study of cosmological 
parameter estimation in the context of extended models 
that encompass a larger model parameter space than the 
standard, flat ACDM cosmology. Using the 6-parameter 
vanilla model as a basis, we include as additional pa- 
rameters only those that are physically motivated, such 
as a nonzero neutrino mass. We consider a 11-parameter 
model and subsets thereof, in contrast with the vanilla-|-l 
approach adopted in most previous analyses which treats 
one extra parameter at a time. 

In this more general framework, we find that in the 
context of standard slow-roll inflation, constraints on the 
dark matter and dark energy parameters can be substan- 
tially altered. If only CMB and LSS data are used, the 
larger parameter space introduces new, strong parame- 
ter degeneracies, e.g., between the physical dark matter 
density flch^ and the dark energy equation of state w. 
These degeneracies can be broken to a large extent by 
adding type la supernova and baryon acoustic oscillation 
data to the analysis. However, even with this expanded 
data set, we find that the bound on the physical dark 
matter density flch^ relaxes by more than a factor of two 
compared to the vanilla model constraint. 

In the same spirit, we have studied how bounds on 
the infiationary parameters n^, r, and as are affected 
by the introduction of extra parameters in the analy- 
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TABLE III: The ID marginalised 95 % C.L. allowed ranges for Qch^ and w for various parameter and data sets. 



Parameter set 


Data set 






■w 




B 


WMAP+SDSS 


0.092 ^ 0.136 


-1.44 




-0.76 


B 


WMAP+SDSS+SNLS+BAO 


0.100 ^ 0.123 


-1.12 




-0.87 


C 


WMAP+SDSS+SNLS+BAO 


0.100 ^ 0.123 


-1.11 




-0.86 


A 


WMAP+SDSS+SNLS+BAO 


0.094 ^ 0.136 


-1.19 




-0.88 



-0.6 
-0.8 
3 -1.0 
-1 .2 
-1.4 

0.04 
0.02 
a 0.00 
-0.02 
-0.04 





sis. We find that the simplest X(p'^ model of inflation is 
still compatible with all present data at the 95 % level, 
in contrast with other recent analyses P, U, H^l- The 
source of this apparent discrepancy is a strong degener- 
acy between the tensor-to-scalar ratio r and the neutrino 
fraction f^, the latter of which was fixed at zero in the 
analyses of P, 0, [13] ■ Reversing the argument, if Xcj)'^ 
is the true model of inflation, then it strongly favours 
a sum of quasi-degenerate neutrino masses between 0.3 
and 0.5 eV, a range compatible with present data from 
laboratory experiments. This represents a clear example 
of how neutrino masses well within laboratory limits can 
bias conclusions about other, seemingly unrelated cosmo- 
logical parameters. 
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